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1 INTRODUCTION 

Observations show that a significant fraction, perhaps a 
majority, of AGN of different types are obscured by a 
screen of a cold dusty matter, thought to be a molec- 
ular torus-like structure with a scale between 1 to a 
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many lines of observational evidence suggest that the un- 
obscured AGN have similar torii as well, but we are viewing 
these AGN at an angle allowing a direct view of their cen- 
tral engines. The orientation-dependent obscuration is thus 
said to unify the differ ent AGN classes iAntonucci fc Miller! 
ll985l:lAntonuccilll993l) . 

Despite the observational importance, the theory of 
molecular torii is still in an explorator y state, which is in- 
dicativ e of the d ifficulty of the problem. [ Krolik fc Beeelmanl 
l|l98Sft (see also lKrolik fc Begelmanlll98a) showed that the 
large geometrical thickness of the torus and yet its ap- 
parently low temperature are consistent only if the torus 
is made of many molecular clouds moving with very high 
random velocities (Mach number £ 30). Given that stellar 
feedback processes appear inefficient to explain these large 
random speeds, the authors were "forced to seek a much 
more speculative solution: viscous heating of the cloud sys- 
tem due to partially elastic cloud-cloud collisions." However, 
until now there is no definitive answer (via numerical simu- 
lations) on whether clouds colliding at large Mach numbers 
would behave even partially elastic, or would rather share 
and dump the random velocity component and collapse to 
a disk configuration. 



ABSTRACT 

Orientation of parsec-scale accretion disks in AGN is likely to be nearly random for 
different black hole feeding episodes. Since AGN accretion disks are unstable to self- 
gravity on parsec scales, star formation in these disks will create young stellar disks, 
similar to those recently discovered in our Galactic Center. The disks blend into the 
quasi-spherical star cluster enveloping the AGN on time scales much longer than a 
likely AGN lifetime. Therefore, the gravitational potential within the radius of the 
black hole influence is at best axi-symmetric rather than spherically symmetric. Here 
we show that as a result, a newly formed accretion disk will be warped. For the simplest 
case of a potential resulting from a thin stellar ring, we calculate the disk precession 
rates, and the time dependent shape. We find that, for a realistic parameter range, 
the disk becomes strongly warped in few hundred orbital times. We suggest that this, 
and possibly other mechanisms of accretion disk warping, have a direct relevance to 
the problem of AGN obscuration, masing warped accretion disks, narrow Fe Ka lines, 
etc. 

IWada fc Normanl faOO'j) took a point with 
iKrolik fc BeeelmarJ il988t) assertion of the weakness 
of the stellar processes. Considering "large", e.g. 100 parsec 
torii, and via numerical simulations of supernova explosions 
resulting from star formation in a self-gravitating, massive 
disk, they have shown that large random cold gas velocities 
result from the interaction of the gas with the supernova 
shells. However as such about 99 % of the gaseous disk in 
the simulations is consumed in the star formation episode 
rather than being accreted by the black hole. Unless most 
of the stars formed in the disk are later somehow accreted 
onto the SMBH, this mechanism can only work for torii on 
scales much larger than the SMBH gravitational sphere of 
influence, Rh, i.e. where the stellar mass is much greater 
than A/bh- It is also not clear whether this mechanism 
would work for smaller accretion/star formation rates be- 
cause then supernova explosions become too rare. Finally, 
time variability of obscuring column depths argues for much 
smaller radial sizes of the "torus", e.g. i n the range 0.1 — 10 
parsec jRisaliti. Elvis, fc Nicastroll2002l) . 

Here we would like to emphasize that accretion disks 
in AGN are very unlikely to be planar. There are several 
mechanisms that are capable of producing strong warps. We 
believe that such warps have to be an integral part of the 
AGN obscuration puzzle. In particular, we shall describe 
and calculate a disk warping mechanism driven by an axi- 
symmetric gravitational potential. 

The sequence of events that takes place in an AGN feed- 
ing cycle in our model is as following. First of all, a merger 
with another galaxy or a satellite, or another source of cold 
gas, fills the inner part of the galaxy with plenty of gas. This 
gas will in general have a significant angular momentum ori- 
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ented practically randomly. The gas will start settling in a 
disk that is to o massive and too cold to be sta ble again st self- 
gravity (e.g., |PaczYnski||l^f<^ |Shlosma n fc BegelmarJll989t 
ICollin fc Zahnlll999l : lGoodmanll2003h . Stars are formed in- 
side the accretion disk, producing a flat stellar system. This 
long-lived axi-symmetric (or perhaps also warped) structure 
will torque any orbit which is not exactly co-planar with it, 
resulting in a precession around the symmetry axis. Now, as 
the time goes on, the orientation of the angular momentum 
vector of the incoming gas changes, and the newly built disk 
is exposed to the torque from the stellar disk remnant. Dif- 
ferent rings of the disk precess at different rates, thus the 
disk becomes warped. 

We shall emphasize that existence of these stellar disk 
remnants is hardly a question based on the severe self- 
gravity problems for the standard accretion disks at large 
radii. The recent observational evidence supports this con- 
clusion too. The best known ex ample of such flat stellar 
systems is in our Galactic Center dGenzel et aljl20o3) where 
two young stellar rings are orbiting the SMBH only a tenth 
of a parsec away. The orbital planes of these rings are 
oriented at very large angles to the Galactic plane. Such 
flat stellar systems are also being found elsewhere (e.g., 
iKrainovic fc Jaffdl 20041) . 

In this paper, we concentrate on the linear gravitational 
warping effect to investigate its main features. Taking the 
simplest case of a warp produced by a massive ring, we cal- 
culate the gravitational warping torque. Starting from a thin 
test accretion disk, we calculate its time dependent shape. 
Under quite realistic assumptions, the disk becomes strongly 
warped in some 10 2 — 10 3 orbital times, i.e. in 10 5 — 10 7 
years. Since the warping is gravitational in nature, gaseous 
and clumpy molecular disks alike are subject to such de- 
formations. We speculate that in the non-linear stage of the 
effect, a clumpy warped disk will form a torus-like structure. 



2 TORQUES IN A LINEAR REGIME 

The basic physics of the effect i s very simple and has been 
explained by e.g.. iBinnevI il992T) in his consideration of the 
galactic warps (most galactic disks are somewhat warped). 
Consider a nearly circular orbit of radius R for a test particle 
in an axi-symmetric potential (e.g. a central point mass plus 
a flat disk). Let the axis z to be perpendicular to the plane 
of the disk. The particle makes radial and vertical oscilla- 
tions with slightly different frequencies, and the difference 
is the precession frequency, u p . Due to the symmetry, the 
z-projection of the angular mom entum of the particle is ex - 
actly conserved (see, e.g., §3.2 in lBinnev fc TremaindlT987t) . 
Thus the angle between the angular momentum of the parti- 
cle and the z-axis, (3, is a constant. The line of the nodes (the 
line over which the two planes intersect), however, precesses. 

Now consider a "test-particle" disk initially in the same 
plane as the circular orbit. The disk is a collection of rings, 
i.e. circular orbits. Since precession rates uj p are different for 
different R, rings turn around the z-axis on unequal angles. 
The initially flat disk will be warped with time. 



2.1 Gravitational torques between two rings 

Let us calculate the gravitational torque between two rings 
with radii Ri and R, inclined at angle (3 with respect to 
each other. We work in two rigid coordinate systems, x, y, z 
and x' , y , z 1 , centers of which are at the SMBH. For conve- 
nience, we place the first ring in the z — plane, whereas 
the second ring is at z = plane. Angle (3 is obviously the 
angle between the axes z and z 1 . Further, the axes x and x' 
are chosen to coincide with the line of the nodes. 

The total torque exerted by the second ring on the first 
one is 



T2i = GaiRiaR 



[fi x f] 

f — fi| 3 



(1) 



where <n = Mi/2wRi and a = M/2nR, with Mi and M 
being the masses of the two rings, respectively. The integra- 
tion goes over angles <j>\ and <f> that are azimuthal angles in 
the respective frames of the rings. From this expression it 
is immediately clear that co-planar rings do not exert any 
torque on each other, [fi x r\ — 0. In addition, if (3 = n/2, 
the integral vanishes as well because for each fi the oppo- 
sites sides of the second ring (e.g. (j) an d <f> + tt) make equal 
but opposing contributions. 

It is also possible to show that due to symmetry the 
torque's z-projection vanish. Thus we only have the T2i, x = 
—t~12,x component, meaning that the angular momentum 
vector of the ring will rotate without changing its magni- 
tude. Also, if Mi S> M, then one can neglect warping of the 
first ring. 

The absolute distance between two ring elements with 
the respective angles (f>i and <j> is 



[fi - rf = R\ + R 2 - 27?ii?cosA , 
where 

cos A = cos (3 sin <f> sin cj>i + cos cj>\ cos <j> 



(2) 



(3) 



Without going into tedious detail, we write the torque ex- 
pression separating out the leading radial dependence and 
the integral over angles, which we label 1(5,(3): 



T21,; 



with 



GM1MR1R 
(R 2 + J? 2 ) 3 / 2 



1(6,(3) , 



1(8,(3) = sin/3 
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3/2 



and 



2fli_R 
Rf + R 2 



(4) 



(5) 



(6) 



Total angular momentum of the second ring is L = 
MQkR 2 - Recall that L z — Lcos(3 — const, whereas the 
component of L in the plane of the first ring (z = 0) pre- 
cesses. We chose it to be initially in the y-direction, so that 
L y (t — 0) = Lsin/3 (cf. equation 1121 below*) . The precession 
frequency for the second ring is then 



Tl2, : 



Tl2,; 



MQ, K R 2 sin (3 ' 



(7) 



In general, the integral in equation |S] is calculated 
numerically, but for 5 1 one can decompose [1 — 
<5cosA]- 3/2 » 1 + 3/25 cos A and obtain 1(5,(3) » 
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(3/8)5 sin 2/3. When this approximation holds, that is when 
R S> Ri or R <C Ri , the precession frequency for the second 
ring is 

3Mi n R?R\ 



0/ 



4M BH 



cos (3 



[R 2 + R\ 



2l5/2 



a maximum at radius R m at which R m /Ri 
maximum growth rate of the warp is thus 

Mi 



max \(jp\ ~ 0.085r2/f cos /3 



Mr 



(8) 



3/7. The 



(9) 



Mi/Mbh may be expected on the order of a percent or so 
since this is when the accretion di sks become self-gravitating 
(e.g., iNavakshin fc CuadralEool . and the resulting stellar 
mass would probably be of that order too. In this case one 
notices that to produce a sufficiently large warp one has to 
wait for at least ~ 10xMbh/Ml ~ lOOOft^ 1 or ~ 150 orbital 
times at the radius of the maximum warp R m . While this is a 
fairly long time, it is still much shorte r than the correspond- 
ing dis k viscous time (cf., e.g., Fig. 2 in lNavakshin fe Cuadral 
J2004M . 

The assimptotic dependence of precession frequency on 
radius R is 

Qk for R <C Ri , and 



3Mi a R* 
- — T~ r — cos 3 — j 

4M BH R\ 



3Mi R\ 

'TT7 cos P T^7 K 

4M BH R 2 



for R > R x 



(10) 



(11) 



A thin massive stellar ring would thus only warp a range of 
radii, leaving the portions of the disk much internal and also 
much external to it unaffected. 



2.2 Test disk warping 

We now want to calculate the shape of a light non self- 
gravitating disk warped by a massive ring Ri. The disk is 
treated as a collection of rings with different radii R and 
negligible mass. We assume that the mass of the whole disk, 
M, is negligible in comparison with the mass of the ring, 
M <C Mi, and that therefore the massive ring's orientation 
(angular momentum) does not change with time. 

In application to the real disks, it should be remembered 
that any attempt to warp a disk will be resist ed by disk vis- 
cous forces (e.g.. Eardeen fc Pettersonlll975f) that transfer 
the angular momentum through the disk. However, we are 
interested in the outermost regions of AGN accretion disks 
where they "connect" to the galaxy. On these, a tenth of a 
parsec to 100 parsec scales (the range depends on M bh), 
the usual o-disk viscosity dShakura fc Sunvaevl I1973T) be- 
comes ineffective. The viscous time scale becomes too long 
and the disks are believed to be prone to self-gravity (e.g., 
| Paczvnskll978l:|Shlosman fc BegelmarJl989T:ICollin fc Zahnl 
ll999l:lGoodmar]l2003h . Therefore we are justified in neglect- 
ing the restoring viscous forces. A much more important 
effect will be the restoring force from the self-gravity of the 
disk that is being warped, but we differ the study of this 
and other non-linear effects to a future paper. 

It is convenient to work with angle /3, already intro- 
duced, and one additional angle, 7. Recall that angle j3 is 
the local angle of the ring's tilt to the z-axis and it remains 
constant in the test particle regime. Angle 7 is needed to 



introduce the projections of the unit tilt vector no rmal to 
the ring, l(R,t), on the x and y axes dPringlelll996l) : 

1= (cos 7 sin j3, sin 7 sin /3, cos 0) . (12) 

Angle 7 therefore describes the precession of each ring 
around the z-axis. With the chosen coordinate system, at 
time i = 0, 7 = 7r/2 (and we get, in accord with conventions 
of L y (0) = Lsin/3, L x (0) = 0). As each of the disk 

rings precesses, 

j = ir/2+ujpt . (13) 

It is useful to write the expression for the unit tilt vec- 
tor in the (x' , y' ,z') coordinate system rigidly bound to the 
initial disk plane: 

l' x = sin j3 cos 7 , 
l' y — — cos (3 sin (3{1 — sin 7) , 



l' z = cos 2 (3 + sin 2 (3 sin 7 



(14) 
(15) 
(16) 

Note that if f3 = 0, l' x = l' y = 0, e.g. the rings are never 
tilted with respect to the z — z' axis, as it should. Also, 
when 7 = 7r/2, I' indeed coincides with the z'-axis, i.e. the 
disk is flat. 

Using equations 1141161 we can now calculate the shape 
of the warped disk in that system given the function y(R, i). 
To accomplish this, one first introduces the azimuthal angle 
(f> on the surface of the ring. The coordinates of the points 
on the ring, r , are then given by equations 2.2 and 2.3 in 
iPringld il996t) . The corresponding coordinates in the primed 
system of reference are easily obtained by x' = (rej'), etc., 
where ej and so on are the unit coordinate vectors of the 
primed system. The result is: 



x 
R 



v_ 

R 



cos 13 cos 7 sin <f> + sin 7 cos <j> , (17) 
cos (3 cos 7 cos 4> + sin </>(cos 2 (3 sin 7 + sin 2 (3) , (18) 
= — sin f3 [cos 7 cos (f> + (1 — sin 7) cos (3 sin <f)] . (19) 



2.3 An example 

To illustrate the results, we calculate the precession rates 
uip(R) for the following case, Ri = 2, Mi = O.OIA/bh, 
(3 = tt/4. The resulting shape of the disk in the original 
un-warped disk plane is plotted in Figure The time is 
in units of fl^ 1 at R = 1, i.e. t = 1 corresponds to time 
t = 1400 r'pi 2 Mg 1 ^ 2 years, where r pc is the distance in units 
of parsec and Ms = Mbh/10 8 Mq. Note that the warp is 
strongest around radius R ~ Ri — 2, as expected. The in- 
ner disk is hardly tilted, which is not surprising given that 
ujp(R) — * for small R (equation 1101 . Same is true for the 
outer radii, where the small tilt of the original plane can be 
seen on the edges of the disk (we do not calculate the tilt 
beyond R = 5.2 and hence the original disk plane is still 
seen on the edges of the Figure). 



3 DISCUSSION 

3.1 Obscuration of the central engine in AGN 

We believe that the gravitational disk warping due to non- 
spherical mass distribution within the SMBH sphere of in- 
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t = 200 




t = 400 




t = 800 




t=1 200 




Figure 1. Snapshots of the shape of a massless accretion disk warped by a stellar ring of radius R = 2 inclined at angle j3 = ir/4 with 

respect to the disk. The snapshots are for four different time as indicated on the top of each panel. The time unit is 1/Qk (R = 1) = 
3/2 1/2 

1400 rp' c M g ' years. At times larger than those used in the Figure, the disk becomes warped so much that, looking from its initially 
non-warped plane, the plane equation z'(x',y') becomes a multiple-valued function for some (x',y'). In reality non-linear effects will 
limit the growth of the warp. 



fiuence is a common occurrence in real AGN. It is hard to see 
why AGN disks should always form in one plane and even 
why they should be planar when they are born iPhinnevI 
1989). The time scales for development of strong warps are 
short or comparable to characteristic AGN lifetimes (which 
are thought to be in the range of 10 7 — 10 8 years). We thus 
expect that the outer edges of the disk will obscure a sig- 
nificant fraction of the sky as seen from the central source 
and h ence be directly r elevant to the unification schemes of 
AGN iAntonuccill99a) . 

Figure [5] shows the warped disk surface from Figure 
at different times as seen from the origin of coordinates. The 
vertical axis shows cos = z'/\/ (x 1 ) 2 + (y') 2 + (z') 2 an d the 
horizontal one shows (j> defined previously. The shape of the 
disk at different times is shown with different colors. At time 
t — the disk is flat and its projection is simply cos 9 — 
for all 4>- At later times the rings of the disk near R = 



Ri = 2 precess and bulge out of the initial disk plane. With 
time the fraction of the sky obscured from the central engine 
becomes greater than a half. In reality non-linear effects and 
interaction of the disk with AGN winds and radiation, etc., 
will become important in shaping the disk for significant 
warps. 



3.2 The non-linear evolution 

In general, the non-linear stage of the evolution of the sys- 
tem of disks or rings of stars and molecular gas is far from 
a thin disk as long as collisions are of minor importance. 
We have already explored the non-linear evolution of col- 
lisionless systems with N-body codes and the results will 
be reported elsewhere. When two massive disks warp each 
other, mixing occurs when 7 — it/ 2 becomes large, and the 
resulting configuration reminds that of a torus. 
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Warped disk as seen from SMBH 




100 200 300 

0, degrees 

Figure 2. The accretion disk surface as seen from the SMBH 
location for four different times: ( = 100, 400, 800, 3200 for the 
black, green, red and yellow dots, respectively. Note that at the 
largest time most of the available solid angle is obscured by the 
strongly warped disk. 



Now, inelastic dissipative collisions between gas clumps 
will eventually become important. Relaxation of a collision- 
less N-body system leads to non-circular orbits, and hence 
different disk rings will start to overlap. Collisions should 
then tend to destroy the random motions and should estab- 
lish a common disk "plane" . However, such a disk will nor- 
mally be warped itself. Therefore inelastic dissipative col- 
lisions do not necessarily turn off the obscuration in our 
model. Further, if molecular gas clumps are constantly sup- 
plied from the outside and come with fluctuating angular 
momentum, the disk may never arrive in a flat thin config- 
ura tion (see also |PjiimieJl989|) . 

iNenkova. Ivezic. fc Elitzurl J2002f) and 

lRls^h!^F^lvrs^^^castrol J2002J) convincingly argued 
that the AGN obscuring medium must not be uniform in 
density. This does not contradict our model at all since 
AGN accretion disks on large (e.g. 0.1 pc and beyond) 
scales are self-gravitating unless th e accretion rates are tiny 
(e.g., IShlosman fc Begelmanl ll989'l. 

The source of warping potential does not have to be 
a thin stellar ring or a disk: it may be any non-spherical 
distribution of stars in the SMBH vicinity that retains a non- 
zero quadrupole moment; it can also be a second (smaller) 
super-massive black hole during a merg er of two galaxies. 

Both the collision-dominated (iKrolik fc Begelmanl 

|l988|) and the stellar-feedback inflated torii 
JWada fc Normanl l2002jl share a common starting point: 
the torus is the result of some internal disk physics. If our 
interpretation applies, AGN torii lend their existence to the 
way in which the cold gas arrives in the central part of the 



galaxy . Compared with the model of IKrolik fc Begelmanl 
(1988), large random speeds for the cold gas are not 
required in our model. Although the gas may be quite 
high up away from the original plane of the disk, the disk 
is locally coherent and thin. High speed elastic collisions 
between molecular clouds are thus not needed to explain 
obscuration of AGN in our model. 

There are also two other mechanisms potentially able 
to produce warped disks at parsec and beyond scales around 
AGN. As mentioned in the X-ray binaries context, accretion 
disks develop twists and warps due to inst abilities driven by 
X-ray heated wind off the disk surface iSchandl fc Meven 
1994). Same can be achie ved by the rad iation pressure force 
from the central source JPringlelll99rl . However, it seems 
that the majority of disks in X-ray binar ies are either not 
warp ed or warped not strongly enough llQgilvie fc Dubusl 
2001) to provide the large obscuration needed for the AGN 
unification schemes. In contrast, the warping mechanism dis- 
cussed here is not applicable to X-ray binary systems, and 
hence it may be natural that AGN disks are stronger warped 
than X-ray binary disks on appropriately scaled distances 
from the center. 



3.3 The stellar disks in Sgr A* 

The two young stellar disks discovered recently in Sgr A* 
present a challenge to the usual star formation modes be- 
cause the gas densities required to avoid tidal shearing are 
many orders of magnitude larger than the highest densities 
observed anywhere in the galactic molecular clouds. On the 
other hand, star formation inside a massive accretion disk is 
a long expected ou tcome of the self-gravitational insta bility 
of such disks (e.g. lPacz vnskilll978t ISh losman fc Begelmanl 
119891: ICollin fc Zahnlll999l: iGoodmanl l2003h . As such, the 
young stars in the Sgr A* star cluster therefore are a first 
example of star formation in this extreme environment. It 
is very likely that the star formation efficiency and the ini- 
tial mass function (IMF) are quite different in the imme- 
diate AGN vicinity and elsewhere in a galaxy. The "astro- 
archeology" of Sgr A* can be used to study these issues. 

The gravitational warping effect discussed in this pa- 
per constrains the time-averaged mass of the outer ring, 
M ou ter, since the moment of its creation, assumed to be 
t = 2 x 10 6 years. The i nner stellar ring is rather well de- 
fined (iGenzel et alJl2003ft . and we thus estimate ui p t for the 
inner ring to be smaller than 7r/4. Taking the radius of the 
inner stellar ring to be R = 3" ~ 4 x 10 4 R s for the GC 
black hole, and for the outer, Ri = 5", and using equation 
13 with cos/3 = 1/4, one obtains M ou tcr < 10 5 M . Prelimi- 
nary numerical N-body simulations show this limit may be 
even smaller. 



3.4 Other implications for AGN 

There are clearly other observational implications of gravi- 
tationally warped disks. For example, Narrow Fe Ka lines, 
observed in many Seyfert galaxies, can be explained with X- 
ray reflection off such warped cold disks. In addition, warped 
disks will yield different coherent paths for maser amplifica- 
tion than flat disks do, which may be a part of the explana- 
tion for the complexity of the observed AGN maser emission. 



6 Nayakshin 



4 CONCLUSIONS 

We argued that clumpy self-gravitating accretion disks in 
AGN are generically strongly warped. We believe such warp- 
ing should be an integral part of the explanation for the 
AGN unification schemes. Our model provides arguably the 
easiest way to obscure the central engine without the need 
to lift cold gas clouds off the disk plane high up via elastic 
collisions, supernova explosions or winds. 
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